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ABSTRACT 

We present an analysis of CO spectroscopy and infrared-to-millimetre dust photometry of 11 bright far-infrared/submillimetre sources discovered 
through a combination of the Planck all-sky survey and follow-up Herschel- SPIRE imaging - “ Planck’s Dusty Gravitationally Enhanced sub- 
Millimetre Sources.” Each source has a spectroscopic redshift z = 2.2-3.6 from a blind redshift search with EMIR at the IRAM 30-m telescope. 
Interferometry obtained at IRAM and the SMA, and optical/near-infrared imaging obtained at the CFHT and the VLT reveal morphologies con¬ 
sistent with strongly gravitationally lensed sources. Additional photometry was obtained with JCMT/SCUBA-2 and IRAM/GISMO at 850//m 
and 2 mm, respectively. The spectral energy distributions of our sources peak near either the 350//m or 500 pm bands of SPIRE. All objects are 
bright, isolated point sources in the 18" beam of SPIRE at 250 pm, with apparent far-infrared luminosities of up to 3 X 10 14 L© (not correcting 
for the lensing effect). Their morphologies and sizes, CO line widths and luminosities, dust temperatures, and far-infrared luminosities provide 
additional empirical evidence that these are strongly gravitationally lensed high-redshift galaxies on the submm sky. We discuss their dust masses 
and temperatures, and use additional WISE 22-pm photometry and template fitting to rule out a significant contribution of AGN heating to the total 
infrared luminosity. Six sources are detected in FIRST at 1.4 GHz. Four have flux densities brighter than expected from the local far-infrared-radio 
correlation, but in the range previously found for high-z submm galaxies, one has a deficit of FIR emission, and 6 are consistent with the local 
correlation, although this includes 3 galaxies with upper limits. The global dust-to-gas ratios and star-formation efficiencies of our sources are pre¬ 
dominantly in the range expected from massive, metal-rich, intense, high-redshift starbursts. An extensive multi-wavelength follow-up programme 
is being carried out to further characterize these sources and the intense star-formation within them. 
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1. Introduction 


The brightest, most strongly gravitationally lensed galaxies in 
the high-redshift Universe have an extraordinary potential to 
advance our understanding of the processes that regulate the 
growth of the most massive galaxies we see today. In particularly 
fortuitous cases, gravitational lensing from massive, intervening 
galaxies or galaxy clusters not only boosts the apparent inte¬ 
grated brightness of high-redshift galaxies by factors up to 20- 
60, but also magnifies their images by similar factors while con¬ 
serving surface brightness). Thereby, they allow us to study the 
fine spatial details of intensely star-forming high-redshift galax¬ 
ies at scales much below lkpc, down to around 100 pc (e.g., 
|Swinbank et al.|[20l0| |2011| ). This is more akin to the scales 
of individual star-forming regions in nearby galaxies than the 
galaxy-wide scales (of order a few kpc) with which we must 
otherwise content ourselves at cosmological distances. 

Since the discovery of th e firs t gravitationally lensed galaxy 
in the optical ( So ucail et al.|1987| the first gravitationally lensed 
quasar had already been discovered by Walsh et al.||1979] ), 
strongly gravitationally lensed galaxies have been identified and 


* Based on observations collected with the Herschel and Planck 
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studied in all wavebands from the ground and in space. Strongly 
gravitationally lensed submm galaxies provide an extraordinary 
possibility to probe individual star-forming regions in the most 
intensely star-forming high-redshift galaxies. They are also very 
promising sources to increase our understanding of how the deep 
gravitational potential well of high-redshift galaxies, their high 
gas fractions and gas-mass surface densities, and feedback from 
star formation and perhaps active galactic nuclei, are setting the 
stage for the intense star-formation at fine spa t ial detail (e.g., 
Danielson et al.||2011t ISwinbank et al.||2011[ |Combes et al.| 
2012 ). 


Massive, dust-enshrouded, and relatively evolved high- 
redshift galaxies are characterized by the bright thermal infrared 
emission from dust heated by intense star formation during their 
rapid primordial growth phase. These galaxies have typical stel- 
lar and dynamical masses of a few times 10 10-11 M© (Smail et al. 


|2004HSwinbank et al.|2006|[Nesvadba et al.|2007| ) and form stars 
at prodigious rates of up to about 1000M o yr _1 , which are un¬ 
parallelled in the nearby Universe. As a population, dusty star- 
burst galaxies may have contributed as much as about half of 
the total ener gy production from star formation at th ese cosmic 
epochs (e.g., |Hauser & Dwek|200H|Dole et al.|2006| . 

Unfortunately, given their importance for our understanding 
of high-redshift star formation and galaxy growth, the large stel¬ 
lar and dynamical masses and short evolutionary timescales of 
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these galaxies make them very rare on the sky. The densities 
of far-infrared and submillimetre selected (FIR/submm) galaxies 
that are bright enough to be good candidates for strong gravita¬ 
tional lensing are only around one every few square degrees for 
sourceswith S s on ~ 100 mJy, with large uncertainties. For exam¬ 


ple, [Lapi et al. (2012) predict about 0.003-0.1 deg -2 for sources 
with S 500 > 400 mJy on the sky, for expected maximal gravi¬ 
tational magnification factors of 20-30 adopted in the models. 
Consequently, we may expect to find only a few of these sources 
on the entire extragalactic sky, in acco rdance with models o f 
the diffuse infrared background light (Bethermin et al. 2 012| ). 
Nonetheless, at luminosities above about 10 iJ L o (corresponding 
to S 500 ~ 100 mJy), they are expected to dominate the integrated 
FIR/submm luminosity function ( [Negrello et al.|2007] ). 

A major breakthrough has been the recent discovery by 
Herschel\^jand the South-Pole Telescope (SPT) of sizeable sets 
of strongly gravitationally lensed submillimetre galaxies, with 
typical magnification factors of a few. These are using the new 
generation of wide-field surveys (several thousands of square de¬ 
grees), probing in particular the range in FIR/submm flux density 
between 100 and 200 mJy (e.g.,|Negrello et al.|2070{|Harris et al. 


2012[ |Vieira et al.||2013; Wardlow et al. 2013; Bussmann et al. 


2013). Most gravitationally lensed sources identified in these 
surveys are magnified by individual massive galaxies at interme¬ 
diate redshifts, producing partial Einstein rings of a few arcsec 
in diameter and magnification factors up to about 10. However, 
the number of sources identified in these surveys above flux 
densities of 200-300 mJy remains rather small. For example, 
[Bussmann et al.[ ( |2013| ) list 30 sources with far-infrared fluxes 
in this range, with 6 sources above S 350 =350 mJy. So far, to our 
knowledge only one source with £350 > 500mJy at 350 yum has 
been published by the SPT collaboration ( Vieir a et al. |2013|) . 
Meanwhile, the brightest source in the Bussmann et al. ([2013 ) 
sample has S 350 =484 mJy. 

The extragalactic Herschel and SPT surveys together cover 
about 4 % of the sky, which highlights the importance of hav¬ 
ing an all-sky survey to identify in a systematic way high- 
redshift FIR/submm galaxies above the 100-300 mJy regime. 
These sources are likely to be amongst the most strongly grav¬ 
itationally lensed galaxies on the sky. Planck is the first all-sky 
survey in the submillimetre with the depth and spatial resolution 
necessary to probe the brightest, and presumably most strongly 
gravitationally lensed high-redshift infrared galaxies observable 
to us. The 90 % completeness limit of the Planck Catalogue of 
Compact Sources (PCCS) corresponds to Lir « 6 x 10 13 L o at 
z = 2 ( [Planck Collaboration XXVIII||2014j >. At these luminosi- 
ties, even all-sky surveys may reveal only very small numbers 
of sources ([Negrello et al.||2010[ |Bethermin et al.||2012[ |Lapi 


|et al.||2012| ). For example, |Herranz et aT7f ( 2013| ) fad only very 

small numbers of high-redshift sources in the Early Release 
Catalogue of Compact Sources (ERCSC) from Planck (. Planck 
Collaboration et al. 2011), primarily blazars; most extragalactic 
sources in the ERCSC are low-redshift galaxies. 

We used photometry derived from Planck- HFI (described 
in [Planck Collaboration I|20T4| |Planck Collaboration VI|2014| 

to identify all compact sources in the Planck maps that have 
colours consistent with being exceptionally bright, dusty, in¬ 
tensely star-forming high-redshift galaxies. We then obtained 
far-infrared photometry of the most promising candidates (the 
“HPASSS” programme [Planck Collaboration Int. XXVII|2014 


1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im¬ 
portant participation from NASA. 


Paper I hereafter), using the Spectral and Photometric Imaging 
Receiver (SPIRE) on board the Herschel space telescope. 
This sample, obtained through Herschel “Must-Do” Director’s 
Discretionary Time, includes 234 of the brightest, rarest sources 
in the sky (one source per several tens of square degrees, Planck 
Collaboration, 2015, in prep.). The sample was deliberately se¬ 
lected to only include sources that do not fall into the large 
Herschel survey fields, and which are new to the literature. 
SPIRE confirms that these sources have the typical submm 
colours of high-z infrared and submm galaxies. Most of the 234 
sources of the HPASSS sample are overdensities of multiple 
galaxies with the typical FIR colours of high-z galaxies, and only 
four are Galactic cirrus clouds. Another small subset are bright 
individual, isolated point sources in the 20" beam of SPIRE, 
consistent with being exceptionally bright, presumably strongly 
gravitationally lensed high-redshift galaxies. An overview of the 
HPASSS sample is given in Paper I. 

Here we present the first results of our multi-wavelength 
follow-up of the 11 brightest of these isolated HPASSS sources 
(“Planck’ s Dusty GEMS”), which can be observed from the 
northern hemisphere. All have flux densities at 350/im mea¬ 
sured with SPIRE that are at least S 350 ~ 300 mJy, well above 
the typical range probed by the SPT and wide-field Herschel 
surveys. The brightest source has S '350 = 1050 mJy. Welikala 
et al. (2015, A&A submitted) discuss another set of weaker 
gravitational lenses observed with Planck and the SPT; how¬ 
ever, the 11 sources we discuss here were previously unknown 
and are substantially brighter. Another source, PLCKERC857 
PLCK_G270.59+58.52, taken from the ERCSC, which also sat¬ 
isfies our selection criteria and fortuitously falls into the H- 
ATLAS survey area, has already bee n discussed by |Fu et al 


(|2012|), |Herranz et al.| (|2013|), and |Bussmann et al.| d2013|) 
HLS 091828.6+514223 at z = 5.2 ( [Combes et al.||2012[|Rawle| 


et al.||2013| ) behind the galaxy cluster Abell 773, has been dis¬ 


covered independently from our_survey as part of the Herschel 
Lensing Survey (Egami et al. [2010] ), but is also included in our 
Planck parent sample (Paper I). It forms a 1.5", near-complete 
Einstein ring with a magnification factor of 9 + 2 (jRawle et al.| 
|20j3). 


Our present analysis has three main goals. Firstly, we de¬ 
termine spectroscopic redshifts and provide empirical evidence 
from millimetre photometry and spectroscopy that these are in¬ 
deed high-redshift sources, which owe their exceptional bright¬ 
ness in the submm to strong gravitational lensing. We also show 
ground-based near-infrared imaging, illustrating that the submm 
sources lie behind overdense regions in the intermediate-redshift 
Universe. Secondly, we characterize their global dust and star- 
formation properties and show that their observed far-infrared 
emission is dominated by star formation, not by powerful ob¬ 
scured quasars. Thirdly, we estimate molecular gas masses from 
their CO line emission, calculate dust-to-gas mass ratios, and 
use their FIR luminosities and molecular gas mass estimates to 
show that they are more akin to very rapidly forming “starburst” 
galaxies, than the more gradually, but still intensely star-forming 
high-redshift galaxies on the “main sequence” (e.g., |Elbaz et al.[ 
2011| ). This is the first in a series of papers about “Planck’s Dusty 
GEMS,” and for this analysis we only use parts of the compre¬ 
hensive data sets we already have in hand. Subsequent publica¬ 
tions will provide lensing model solutions, discuss the spatially- 
resolved properties of our sources in a number of wavebands, 
and use the multiple millimetre line detections to investigate the 
detailed gas and star-formation properties of these systems. 

The outline of the paper is as follows. After the descriptions 
of our sample selection in Sect. [2j and our follow-up observa- 
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Fig. 1 . Near-infrared J- band images of the nine sources for which we already have optical/NIR imaging and millimetre interferom¬ 
etry (Fig. [2]). Contours show the Herschel- SPIRE morphology at 250 jam. 


tions and data reduction in Sect. [3j we provide the redshifts in 


Sect. 4.2 and the FIR luminosities, dust masses and temper¬ 


atures in Sect. 0 In Sect. [6] we present the evidence that our 
sources are very strongly gravitationally lensed galaxies, includ¬ 
ing arcsec-resolution submm and millimetre interferometry of 
the dust emission, as well as NIR/optical imaging. In Sect. [7] 
we impose additional constraints from the Wide-field Infrared 
Survey Explorer, WISE and the VLA FIRST survey to demon¬ 
strate that the far-infrared spectral energy distributions of our 
sources are not dominated by radiation from powerful AGN, be¬ 
fore turning to their gas and star-formation properties in Se ct. |9| 
and [10| respectively. We summarize our results in Sect. 
Throughout the paper we adopt a flat Ho = 70 km s -1 Mpc -1 
cordance cosmology with Dm = 0.3 and Da = 0.7. 


11 


con- 


2. Sample selection 

In Fig. [I] we show 9 of our 11 newly discovered gravitational 
lens candidates, which are the brightest amongst our sources se¬ 
lected from the nominal data release of the Planck all-sky sur¬ 
vey. The parent sample consists of two subsets that were iden¬ 


tified with similar colour cuts in the 857-GHz to 545-GHz flux 
density ratio of S ^ 57 /S 545 < 1.5-2. This is well-matched to the 
expected spectral shape of the far-infrared continuum of dusty 
starburst galaxies at redshifts z ^ 2. Six of our sources were 
taken from the PCCS, which includes all sources with S/N > 4 at 
545 GHz on the cleanest 52 % of the sky ( [Planck Collaboration 
|XXVIII|2014| ). Five sources come from a dedicated, blind search 
for high-redshift candidates in the Planck maps, which probes 
fainter sources in the cleanest 35 % of the sky after subtracting 
estimates of the cosmic microwave background and Galactic cir¬ 
rus emission. Table [I] lists the origin of each target. The second 
subsample will be described in detail by Planck Collaboration 
2015 (in prep.). A comprehensive summary of the selection and 
in particular the cleaning algorithm adopted to identify this sub¬ 
sample is given in Paper I. 

All sources were followed up with Herschel- SPIRE photom¬ 
etry at 250 pm, 350 pm, and 500 pm as part of the HPASSS sur¬ 
vey, mostly during “Must-Do” Director’s Discretionary Time, 
and is a subset of the sample presented in Paper I. SPIRE has 
about 10 times greater depth and 20 times higher spatial resolu¬ 
tion than Planck, with beam FWHMs of 18", 24", and 35" at 
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Fig. 2. Near-infrared J -band images of the nine sources from Fig. [I] Contours show the SMA morphology at 850yum. The FWHM 
beam size of all SMA data is about 2" (the FWHM beam size is shown in the bottom left comer of each panel), except for 
PLCK_G080.2+49.8, where it is 0.8". 


250/rm, 350yum, and 500//m, respectively Therefore, SPIRE 
can be used to identify and localize the counterparts of Planck 
sources in the same wavelength regime. We identified all po¬ 
tential gravitationally lensed galaxies as bright, isolated point 
sources in the SPIRE images, with spectral energy distributions 
peaking at flux densities above 300 mJy in either the 350 yum or 
the 500 yum band as observed with SPIRE (in the final calibration 
of our SPIRE maps, one source has S 350=294 mJy). This is well 
matched to the 90 %_completeness_ limit of t he PCCS of about 
600 mJy ( [Planck Collaboration XXVIII|2014| ), and purposefully 
excludes similar sources already identified in other surveys. 

We find a total of 15 gravitational lens candidates. For the 
purposes of t his paper, we exclude the two already described in 
the literature ( |Fu et al.||2012[ |Combes et al.||2012[ see Sect. [T), 
and anot her two that are in the far South. The source described 
by Combes et al. ( 2012| ) is at a higher redshift (z = 5.2) than the 
sources we discuss here. Apart from that, including those other 
sources would not change our general conclusions. 


2 See the SPIRE Handbook 2014, http: //herschel. esac. esa. 
int/Docs/SPIRE/spire_handbook.pdf 


In the present analysis we discuss the 11 Northern sources 
that are already spectroscopically confirmed to be at high red- 
shifts. It is worth noting that by focusing on single, very bright 
sources in the Herschel images, we may have missed sets of 
multiple, but fainter gravitationally lensed objects behind the 
same intervening structures; these we would have identified as 
overdensities of high-redshift infrared galaxies in the overall 
HPASSS sample (Paper I). Irrespective of this, given the extraor¬ 
dinary sky coverage of our parent sample and the brightness of 
our targets, six of which reach or even exceed the completeness 
limits of the PCCS in the 353, 545, and 857 GHz bands of 330, 
570, and 680 mJy, respectively (Table 1 of 
XXVIII||2014| ), the 11 sources discussed 
amongst the brightest individual high-redshift galaxies on the 
sky. 


Planck Collaboration 
here are likely to be 
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Table 1. Mid-infrared to millimetre photometry. The columns are: source name; indication of whether the sources come from the 
PCCS or our own high-z selection (OT2, see Sect. [2]); RA and Dec coordinates; redshift of the high-z source; SDSS redshift of 
the intervening foreground source; WISE 22-pm flux density; Herschel- SPIRE 250-pm flux density; Herschel- SPIRE 350-/im flux 
density; Herschel- SPIRE 500-/im flux density; JCMT/SCUBA-2 850-/mi flux density; and IRAM 30-m/GISMO 2-mm flux density. 


Source 

Survey 

RA 

(J2000) 

Dec 

(J2000) 

^source 

£fg,sdss 

s 22 

[mJy] 

s 250 

[mJy] 


s 350 

[mJy] 


s 500 
[mJy] 

^850 

[mJy] 

S 2000 

[mJy] 

PLCK_G045.1+61.1. . . 

OT2 

15:02:36.04 

+29:20:51 

3.4 

0.56 a ’ b 

<5.4 

161+ 

5 

328+ 

6 

397+ 

6 

130+13 

13 

±2 

PLCK_G080.2+49.8 . . . 

OT2 

15:44:32.40 

+50:23:46 

2.6 

0.5 

<5.4 

220+ 

9 

340+ 

6 

314+ 

6 

87+ 5 

7 

±2 

PLCK_G092.5+42.9 . . . 

PCCS 

16:09:17.76 

+60:45:21 

3.3 

0.5 

<5.4 

765+ 

7 

865+ 

8 

696+ 

7 

191+16 

8.4+0.6 

PLCK_G102.1+53.6. . . 

OT2 

14:29:17.98 

+59:21:09 

2.9 

0.74 a 

8.1+0.7 

327+ 

4 

410+ 

5 

339+ 

5 


5.5+0.7 

PLCK_G113.7+61.0. . . 

PCCS 

13:23:02.88 

+55:36:01 

2.4 

0.3 

6.0+0.7 

672+ 

7 

603+ 

6 

373+ 

5 

108+10 

4 

±1 

PLCK_G138.6+62.0 . . . 

PCCS 

12:02:07.68 

+53:34:40 

2.4 

0.4 

6.2+0.8 

619+ 

6 

664+ 

8 

474+ 

6 

123+14 

7 

±1 

PLCK_G145.2+50.9. . . 

PCCS 

10:53:22.56 

+60:51:49 

3.6 


<5.4 

453+ 

5 

719+ 

7 

781 + 

8 

360+23 

38 

±4 

PLCK_G165.7+67.0 . . . 

PCCS 

11:27:14.60 

+42:28:25 

2.2 

0.34 a ’ b 

10.4+0.9 

867+ 

8 

753+ 

6 

472+ 

5 

90+ 9 

8 

±1 

PLCK_G200.6+46.1. . . 

OT2 

09:32:23.67 

+27:25:00 

3.0 

0.6 

<5.4 

209+ 

4 

294+ 

4 

273+ 

5 

110+10 

7 

±2 

PLCK_G231.3+72.2. . . 

OT2 

11:39:21.60 

+20:24:53 

2.9 

0.1 

<5.4 

299+ 

4 

401 + 

5 

341 + 

6 

111+12 

9 

±1 

PLCK_G244.8+54.9. . . 

PCCS 

10:53:53.04 

+05:56:21 

3.0 

0.13 

<5.4 

1050+10 

1054+10 

777+ 

7 

198+11 

19 

±2 

HLS-J0918 0 . 


09:18:28.6 

+51:42:23 

5.2 



85+ 

8 

168+ 

8 

203+ 

9 

125+ 8 

15 

±7 

HATLAS-J1146 d . 


11:46:37.9 

-00:11:32 

3.3 



323+24 

378+28 

298+24 

93+12 

38 

±6 


a Spectroscopic redshift of intervening foreground source from the SDSS. 

b Multiple foreground sources at a common spectroscopic redshift near the line of sight. 

c This s ource has been found with the Herschel Lensing Survey and has previously been discussed by 


Combes et al. 

|2013| ). The last column lists the 1.2-mm flux density measured with MAMBO on the IRAM 30-m telescope. 


( 2012} and 


Rawle et al. 


This source from the ERCSC (Planck Collaboration et al 
by 


Fu et al. 


2011 


falls serendipitously into the HATLAS field and has been previously discussed 


(2012). The last column lists the 1.2-mm liux density measured with MAMBO on the IRAM 30-m telescope. 


3. Photometry 

3.1. Herschel-SPIRE FIR photometry 


We base our analysis on the Herschel- SPIRE photometry at 
250/im, 350/im, and 500/mi, as well as on ground-based 
JCMT/SCUBA-2 850-yum and IRAM/GISMO 2-mm single-dish 
photometry. For sources in the range z = 2.2-3.6, this corre¬ 
sponds to a wavelength range 50-80 pm and 400-600 pm in the 
rest-frame, which cover the expected peak and Rayleigh-Jeans 
tail, respectively, of the warm dust emission heated by intense 
star formation. 

Herschel-SPIRE observations were carried out between 
December 2012 and March 2013. The Herschel-SPIRE photom¬ 
etry was obtained using STARFINDER ( [Diolaiti et al.|2000| as 
part of the HPASSS survey, and is discussed in detail in Paper I. 
STARFINDER integrates over the point spread function obtained 
directly from the image, and is therefore more suitable for data 
dominated by confusion noise than classical aperture photome¬ 
try. Measured flux densities are between 294 mJy and 1054mJy 
at 350/mi and between 270mJy and 800 mJy at 500 pm, re¬ 
spectively (Table [I]). Typical uncertainties are of the order of 
5-10 mJy, which is close to the confusion level in the maps 
( [Nguyen et al.|2010[ ), but only includes the measurement error, 
not the systematic uncertainty of 7 % inherent in the SPIRE pho¬ 
tometry. However, this has no impact on our results. 


3.2. JCMT/SCUBA-2 photometry at 850pm 

The SCUBA-2 ( [Holland et al.||2013[ ) data were taken between 
September 2012 and May 2014 in moderate conditions with in¬ 
dividual observing times of 15 min per source. The data were re¬ 
duced using a configuration file optimized for point-source cal¬ 
ibrators using the smurf data reduction software package for 
SCUBA-2 ( [Chapin et al.|[2013| ). Flux densities were extracted 


using aperture photometry with a 30" diameter aperture, where 
the background was estimated within an annulus with inner and 
outer diameters of 37.5" and 60", respectively. Uncertainties are 
between 4 mJy and 21 mJy per beam. These flux densities are 
then corrected for missing flux density due to the aperture size 
by dividing by 0.85, as described by Demps ey et al.| ([2013 ). An 
8 % calibration uncertainty is added in quadrature to the photo¬ 
metric errors. The resulting 850/tm flux densities are listed in 
Tabled] 


3.3. SMA 850pm interferometry 

All sources were also observed in the continuum with the SMA 
in the 850 yum band. Observations were carried out between 
June 2013 and June 2014. Data were taken under good to ex¬ 
cellent conditions with pwv<2 mm in shared-track mode to ob¬ 
tain good uv coverage, in spite of observing each source with 
less than one track; this was made possible by their extraordi¬ 
nary brightness in the submm. Integration times per source are 
between 2 and 7 hrs. All sources but PLCK_G080.2+49.8 were 
observed through programme 2013B-S050 in the compact con¬ 
figuration, with a beam of about 2" x2". PLCK_G080.2+49.8 
was observed through DDT programme 2013A-S075 in the ex¬ 
tended configuration, giving a beam of 0.8" x0.5". Data were 
calibrated in IDL using the MIR package, and analysis and imag¬ 
ing utilized the MIRIAD package. A full discussion of the inter¬ 
ferometry is beyond the scope of the present analysis, and will 
be presented in a subsequent publication. Here we only use the 
850 pm morphologies to illustrate that these are indeed strongly 
gravitationally lensed galaxies. Comparison with the single-dish 
flux densities from SCUBA-2 suggests that we recover at least 
80-90 % of the total flux density, implying that we have typically 
not missed fainter, more extended components. 
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3.4. I RAM-30-m /GISMO 2-mm photometry 

Two-mm continuum observations were carried out with the 30- 
m telescope of IRAM between 17 and 23 April with 1.3-1.7 mm 
of precipitable water vapor, and between 29 October and 5 
November, 2013, with 3.4-8.6 mm of preciptable water vapor 
(programmes 222-12 and 100-13). 

We used the 8x16 pixel bolometer camera GISMO, which 
covers a frequency range of 140-162 GHz ( [Staguhn et al.|2012| 
in the 2-mm band. At the 30-m telescope of IRAM, GISMO cov¬ 
ers a 1.8' x 3.7' field of view with a 21" x 21" beam. We used 
2' x 2' Lissajous maps with a relative flux density stability of 
about 8 %, which is optimized to obtain high signal-to-noise ra¬ 
tios of relatively faint objects like ours, with only a small number 
of bad channels. Total integration times per source were between 
10 and 100 minutes, depending on the expected flux density of 
the target. 

We used the CRUSH software package (Kova cs]|2013[ ) to re¬ 
duce and calibrate individual scans, and to combine them into the 
final image. We used the “faint” option, which is well adapted to 
signals around lOmJy, detected at signal-to-noise ratios of less 
than 10 per scan. We detected all targets as point sources, with 
S/N between 3.4 and 14.8. 

3.5. CFHT and VLT optical/NIR photometry 

In order to identify and broadly characterize the foreground lens- 
ing structure, we also obtained near-infrared optical imaging 
of our sample using Megacam and WIRCAM on the Canada- 
France-Hawaii Telescope (CFHT), and with HAWKI and FORS 
on the Very Large Telescope (VLT) of the European Southern 
Observatory (ESO). Here we use the optical imaging to high¬ 
light that our sources lie behind massive intervening structures, 
either galaxy groups or clusters. Observations of the overall sam¬ 
ple are still on-going, but have already provided us with optical 
or NIR imaging of all but one target in at least one band. We 
used the scamp and swarp software (Bertin [2010a|b| ) to register 
our images relative to the USNO-B2 catalog, with a typical po¬ 
sitional uncertainty of 0.2"-0.5". A full discussion of these data 
and what they imply more quantitatively for the lensing config¬ 
uration will be presented in a subsequent publication. 


4. Blind spectroscopic redshift survey in the 
millimetre 

4.1. IRAM 30m/EMIR spectroscopy 

We performed a blind redshift search in the 3-mm and 2-mm 
bands for all 11 targets using the wide-band heterodyne receiver 
EMIR at the 30-m telescope of IRAM. Following a pilot pro¬ 
gramme to measure a spectroscopic redshift of our first source, 
PLCK_G080.2+49.8, through regular programme 82-12 and 
Director’s Discretionary programme D05-12, we obtained an¬ 
other 75 hrs of observing time through Director’s Discretionary 
programme D09-12 and the regular programme 094-13 in April 
and June 2013. For all sources we used the WILMA and FTS 
backends during good to variable conditions. 

Individual scans were 30 seconds long, and we observed sets 
of 12 scans followed by a calibration. Data were reduced us¬ 
ing CLASS ( [Gildas Team|2013| ). We took advantage of a dedi¬ 
cated routine kindly provided by C. Kramer to individually cor¬ 
rect the baselines in each of the sub-bands of the FTS backend. 
The full set of lines will be discussed in a subsequent publica¬ 
tion; here we only use the lowest-/ CO transition available in 


each source, which provides additional empirical evidence that 
these are strongly gravitationally lensed galaxies. The redshifts 
are also derived from these lines, and found to be consistent with 
the full set of available lines per galaxy. 


4.2. Spectroscopic redshifts 


To obtain a spectroscopic redshift for each source, we started 
with a blind line search in the 3-mm atmospheric window, which 
we cover almost entirely with two interleaved tunings of EMIR 
centred at 89.4 GHz and 97.4 GHz, respectively. We used the 
WILMA and FTS backends in parallel, which have band widths 
of about 4 GHz and 8 GHz, respectively. 

We discovered a bright emission line in each source in one 
of the tunings in the 3-mm band (Fig. [3]). Subsequently we cal¬ 
culated all possible redshifts compatible with the observed fre¬ 
quency of the line, and tested our redshift hypothesis by search¬ 
ing for a second line (typically CO(5-4) or CO(4-3)) at the pre¬ 
dicted higher frequency. This required a separate frequency tun¬ 
ing, with the choice depending on each redshift hypothesis (typ¬ 
ically lying in the 2-mm band). We started with the redshift that 
was closest to the photometric redshift derived from our SPIRE 
photometry, assuming a dust temperature of 7^ = 30 K, which 
is the temperature of the dust component that dominates the 
FIR spectral energy distribution of the Cosmic Eyelash (Ivison 
|et al.|2010bj ). This yielded the correct redshift in all but the two 
galaxies, specifically those with the highest dust temperatures 
(Table [2]). In the end, our EMIR follow-up spectroscopy led to 
accurate and secure spectroscopic redshifts for all 11 targets, 
with 2-8 lines detected per source; this provides a wide range of 
constraints on the physical properties of the gas in these galaxies, 
which we will discuss in more detail in forthcoming papers. 


Planck's Dusty GEMS fall in the redshift range of z = 2.2- 
3.6, comparable to that of radio-selected submillimetre galaxies 
in the field ( [Chapman et al.||2005] ), and very similar to that of 
gravitationally lensed submm galaxies in the H-ATLAS survey 
(z = 2.1 -3.5; [Harris et al.|2012[|Bus smann et al.|2013"] l. The av¬ 
erage redshift of Planck's Dusty GEMS, z — 2.9 ± 0.4 (we give 
the width of the distribution here), is somewhat lower than the 
redshift range of the bright gravitationally lensed submm galax¬ 
ies from the SPT survey. The SPT sample has a mean redshift of 
z - 3.5, obtained from at least two lines for 12 of their 26 targets, 
and a combination of single-line detecti ons and FIR p hotomet¬ 
ric constr aints for the remaining targets ( Vieir a et al.|2013| WeiB 


et al.|2013] ). One reason that we probe somewhat lower redshifts 


may be that we select our sources at shorter wavelengths (350- 
850yum, compared to 1.4 and 2.0 mm for the SPT). However, 
we stress that our parent sample from Planck does include 
sources with higher redshifts, such as HLS 091828.6+514223 
at z = 5.2 ( [Combes et al.|2012l |Rawle et al.|2013| ), a source we 
selected independently from its confirmation as a high-z source 
through the HLS survey. We also note that the Planck selec¬ 
tion is based on 350/rm, 550/rni, and 850/rm measurements, 
covering somewhat longer wavelengths than the blind Herschel 
surveys. The Planck high-z sample has a range of FIR colours, 
and the redshift distribution derived from the subsample of very 
bright gravitational lenses does not necessarily correspond to 
the redshift distribution of the overall Planck high-z sample 
(HLS 091828.6+514223 at z = 5.2 being a case in point). 
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Fig. 3. Lowest-/ CO emission line in each of our targets. CO(3-2) and CO(4-3) are indicated in yellow and green, respectively. Red 
dotted lines show simple single-Gaussian fits. 

Table 2. Fitted parameters to the FIR dust continuum and dust properties. The columns are: best-fit wavelength of the peak of 
the blackbody in the FIR; peak flux density at T max of the blackbody; FIR luminosity integrated from 8-1000yum, including the 
magnification factor p from our modified blackbody fit, and neglecting the possible presence of additional, hotter dust components; 
FIR luminosity derived from DecompIR, and including the fiducial mid-infrared flux from star formation (see Sect.[7]); star-formation 
rate, including magnification factor yu; dust temperature; dust mass, including the magnification factor ji\ and AGN contribution to 
the FIR luminosity as obtained with DecompIR, where such a component was found. 


Source 

^max 

[yum] 

S max 

[ mJy] 

^Lfir 

[10 13 l g ] 

uL tpl 

[10 13 L o ] 

yUSFR 

[M 0 yr -1 ] 

T d 

[K] 

1uM d 

[10 9 m o ] 

^-FIR.AGN / ipIR.tot 

[%] 

PLCK_G045.1+61.1 _ 

. . 420+10 

360+10 

8.4+0.1 

12.1+0.1 

14462+172 

36.0+0.5 

6.5+0.1 


PLCK_G080.2+49.8_ 

. . 370+10 

345+ 5 

4.6+0.1 

7.1+0.1 

7920+172 

33.0+0.5 

7.2+0.2 


PLCK_G092.5+42.9_ 

. . 322+ 5 

900+10 

24.8+0.2 

34.8+0.2 

42698+344 

50.1+0.4 

5.3+0.1 

10 

PLCK_G102.1+53.6_ 

. . 333+ 4 

410+ 5 

7.9+0.1 

11.9+0.1 

13601+172 

41.1+0.3 

4.3+0.1 


PLCKLG113.7+61.0_ 

. . 255+20 

700+30 

9.9+0.2 

12.6+0.1 

17044+344 

45.0+0.4 

3.5+0.1 


PLCK_G138.6+62.0_ 

. . 305+10 

690+ 5 

9.0+0.1 

13.5+0.1 

15495+172 

38.7+0.3 

6.4+0.1 


PLCK_G145.2+50.9_ 

. . 400+10 

785+30 

21.8+0.2 

30.1+0.2 

37533+344 

40.5+0.4 

11.0+0.1 

33 

PLCK_G165.7+67.0_ 

. . 265+ 2 

875+ 3 

10.3+0.1 

13.4+0.1 

17733+171 

42.5+0.3 

5.1+0.1 


PLCK_G200.6+46.1 _ 

. . 350+ 3 

295+ 2 

5.7+0.1 

8.2+0.1 

9813+172 

37.5+0.5 

4.3+0.1 

22 

PLCK_G231.3+72.2_ 

. . 350+10 

402+ 5 

7.5+0.1 

10.7+0.1 

12913 +172 

39.3+0.4 

4.7+0.1 


PLCK_G244.8+54.9_ 

. . 300+ 2 

1135+ 2 

26.5+0.2 

36.3+0.2 

45625+344 

50.0+0.4 

5.7+0.1 
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4.3. Line profiles and luminosities 

Fig. [3] shows the lines we detected in the 3-mm band. Many 
sources exhibit complex line profiles, which can also be seen 
in the higher frequency data (Canameras et al. 2015b, in prep.). 
These may either originate from several, gravitationally lensed 
regions in the same galaxy that are blended in the large (20") 
beam of the IRAM 30-m telescope, or else they represent sets 
of nearby, perhaps interacting galaxies. Alternatively, they may 
represent intrinsically complex gas kinematics in single star¬ 


forming environments, driven by the interplay of galaxy rotation, 
feedback, or perhaps turbulent motion within star-forming re¬ 
gions. High-resolution follow-up interferometry is currently be¬ 
ing analysed to further elucidate their nature. 

We fitted the emission lines shown in Fig. [3] with sin¬ 
gle Gaussian profiles to obtain the line FWHMs and inte¬ 
grated fluxes listed in Table [3] Integrated fluxes are Iqo = 
(7-34) Jy kms -1 , for FWHM = 213-685 km s" 1 . A detailed 
analysis of the spectral properties of our sources, using our 
full sets of spectral line observations, and taking into account 
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changes in line ratios in individual, kinematically distinct emis¬ 
sion line regions, will be presented elsewhere; here we only 
use the spectroscopic redshifts, and the fluxes and FWHM line 
widths of the lowest-/ CO transition observed in each source in 
order to demonstrate that these are indeed gravitationally lensed 
galaxies with extreme magnification factors. 

In Table [3] we also compile luminosity estimates for the CO 
lines that we detect. We follow [Solomon et aL] ( |1997| ) in translat¬ 
ing the integrated line fluxes into CO luminosities (in brightness 
temperature units), setting 

L' = 3.25x10 7 / co (y^) 2 D 2 l (1+z)-\ (1) 

Here the measured frequency v is given in GHz, z is the redshift, 
£>l is the luminosity distance, the CO luminosity L' is measured 
in Kkms" 1 pc 2 , and the integrated CO line flux is in Jykms -1 . 
In Sect. [9] we will also derive molecular gas mass estimates and 
compare with the dust masses obtained from our SED fitting. 

5. Dust properties 

To further characterize our sources, we fitted their FIR-to- 
millimetre photometry with modified blackbody distributions 
between 250 /am and 2000 jam. using the Python curve_fit rou¬ 
tine from the scipy packag^] The results are shown in Fig. [ 4 ] 
Most photometry points were observed with roughly similar 
beam sizes between 15" and 30", and we only see single, very 
bright components in each image, which makes us confident that 
uncertainties related to confusion and multiple sources within 
the same beam do not dominate our photometry. We do see the 
foreground lensing sources at short wavelengths, between the 
optical and near-infrared, including the blue channels of WISE. 
However, their relatively blue colours and locations suggest that 
they do not contribute significantly to the long-wavelength emis¬ 
sion. 

In the restframe of our targets, our data cover the peak of 
the SED and extend well into the Rayleigh-Jeans tail of the 
dust emission, which is a particularly clean probe of the thermal 
emission of the coldest dust component, dominating the over¬ 
all mass budget. We followed, e.g., lBlain et al.|p003| to allow 
for a frequency-dependent optical depth of the FIR emission, 
which we model in the standard way through a power law with 
emissivity index =2.0 and a critical wavelength of 100yum, at 
which the dust opacity becomes unity. This specific choice of 
is somewhat arbitrary, but at the signal-to-noise ratio and sparse 
frequency sampling of our SEDs, the precise value has a minor 
impact (for plausible values of )3 between 1.5 and 2.0) for our 
study. For example, we find typical offsets in dust temperature 
A7d « 2 K between fi = 1.5 and ft = 2.0, comparable to our other 
uncertainties. 

These fits correspond to dust temperatures of T& = 33-50 K, 
with no significant trend with redshift (Table [2]). Uncertainties 
were derived from Monte Carlo simulations, where we varied 
the measured flux density in each band 1000 times within the 
measurement uncertainties, before fitting the SED in the same 
way as for the data. 

The high infrared luminosities and dust temperatures of 
Planck's Dusty GEMS suggest that most of the radiation in our 
sources comes from dense star-forming clouds, which is why we 
consider the optically thick case more relevant to our analysis. 
We note, however, that had we fit for optically thin emission, we 

3 http://docs.scipy.org/doc/scipy/reference/generated/ 
scipy.optimize.curve_fit.html 


would have derived temperatures between T& =30 K and 41 K, 
lower by 3 K to 11 K in each individual source. Fig.[4]shows that 
all sources are adequately fit with a single temperature compo¬ 
nent, in spite of a slight, systematic underestimate of the 850 /i m, 
or alternatively, slight overestimate of the 500-jam flux density. 

The temperature range we find does not differ significantly 
from those of other strongly lensed samples of high-redshift 
dusty star-forming galaxies selected in the FIR and (sub-)mm. 
This is illustrated in Fig. [5] It is interesting that the temperature 
range of our sources includes those of typical starburst galaxies 
at high redshift, as well as the lower temperatures of more grad¬ 
ually, but still intensely, star-forming galaxies along the main se¬ 
quence (e.g., |Elbaz et al.|2011| ). We do not find an obvious trend 
of dust temperature with redshift. This may either be due to the 
small number of targets observed over a relatively small redshift 
range, or because we are sampling environments with a range of 
properties within individual galaxies. 

We derive infrared luminosities by integrating over our best- 
fit single-component modified blackbody SEDs between 8 /am 
and 1000/rm in the rest-frame, using the flux densities given in 
Table [I] and dust temperatures in Table [2] Luminosities thus ob¬ 
tained are in the range Lfir = (4.6-26.5) x 10 13 yu -1 L o , includ¬ 
ing the unknown gravitational magnification factor fi. Following 
Kennicutt (1998 ) we set SFR[M 0 yr 1 ] = 4.5 x KT 37 L F ir [W], 
finding star-formation rates SFR = 8000^16000 [T 1 M©yr _1 , 
suggesting high intrinsic star-formation rates and high magni¬ 
fication factors. Corresponding apparent dust masses are M& = 
(3.5-ll)x 10 9 /r _1 M Q (Table |5J). 

Fig.[4]shows that we cover the peak of the dust SED for most 
sources, but nonetheless our sensitivity to warmer dust compo¬ 
nents probed in the FIR shortwards of 250 jam (60-80 /am in 
the rest-frame) is limited, and we are likely to miss emission 
from warmer dust components. Through template fitting with 
DecompIR (Mullaney et al. 2011| , which we will further dis¬ 
cuss in Sect. [7] we derived fiducial correction factors between the 
modified blackbody fits obtained here, and the templates we used 
to accommodate our FIR-to-submillimetre photometry (and con¬ 
straints from WISE at 22 jam within the same SED fit). We thus 
find that our modified blackbody fits may underestimate the IR 
luminosity of our targets by 10-30%. We also list the respec¬ 
tive values obtained from DecompIR in Table [2 Taken together, 
the direct estimate from the modified blackbody fits, and those 
from the templates, set a plausible range of luminosities. These 
luminosities would need to be corrected towards larger values if 
an additional warmer dust component was found with infrared 
constraints shortward of 250 jam. 

We follow, e.g., Greve et al. { 2012} in estimating dust 
masses, Mj, by setting 


D 2 S 

M d =n~ l 7T f 7° (B Vr (T d ) - B Vr (TcMB(z )))- 1 

(1 +Z)Ky r 


( 2 ) 


where Dl is the luminosity distance to redshift z, S vo the flux 
density at the observed frequency v r (l +z) -1 , which we set to 
600 GHz, (corresponding to the 500 /am band of SPIRE). Tcmb 
is the CMB temperature at redshift z. We used k v J [m 2 kg ~^ = 

0.045 x (v r / 250 GHzf , with p = 2.0 as already stated above. 
With this approach and these assumptions, we find dust masses 
of Md = (3.5 - 11) x 10V -1 M o , with the unknown gravitational 
magnification factor yu. Results for individual sources are listed 
in Tabled 
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Table 3. CO line properties obtained with EMIR on the IRAM 30-m telescope. The columns are source name; CO transition; 
integrated line flux; line luminosity \i L' (in brightness temperature units) derived using Eq.[lJ molecular gas mass (see Sect. [9] for 
details); and Full width at half maximum of the emission lines. 


Source 

Transition 

Line flux fi Ico 

Line luminosity yu L' 

P 3^mol 

FWHM 



[Jykms- 1 ] 

[10 11 Kkms -1 pc 2 ] 

[10 11 M 0 ] [km s' 

- 1 ] 

PLCK_G045.1+61.1 _ 

4-3 

22.9+0.3 

6.9+0.1 

5.6+0.1 

213+ 

11 

PLCK_G080.2+49.8_ 

3-2 

9.2+0.5 

3.2+0.2 

2.6+0.2 

265+ 

12 

PLCK_G092.5+42.9_ 

4-3 

34.3+0.2 

9.7+0.1 

7.8+0.2 

453+ 

3 

PLCK_G102.1+53.6_ 

3-2 

5.7+1.8 

2.4+0.7 

2.0+0.2 

252+ 

10 

PLCK.G113.7+61.0_ 

3-2 

16.5+0.2 

5.0+0.1 

4.0+0.2 

528+ 

5 

PLCK_G138.6+62.0_ 

3-2 

22.0+1.1 

6.8+0.3 

5.4+0.1 

514+ 

40 

PLCK_G145.2+50.9_ 

3-2 

21.9+0.8 

12.7+0.6 

10.2+0.2 

685+ 

17 

PLCK_G165.7+67.0_ 

3-2 

25.4+0.3 

6.8+0.1 

5.4+0.1 

576+ 

4 

PLCK_G200.6+46.1 _ 

3-2 

11.2+0.1 

4.9+0.1 

3.8+0.1 

458+ 

9 

PLCK.G231.3+72.2_ 

3-2 

9.4+0.2 

3.9+0.1 

3.0+0.1 

257+ 

8 

PLCK.G244.8+54.9_ 

3-2 

7.4+1.0 

3.3+0.4 

2.6+0.2 

382+ 

10 



Fig. 4. Modified blackbody fits to the FIR-to-millimetre photometry of our sources as obtained with Herschel -SPIRE at 250, 350, 
and 500yum, with JCMT/SCUBA-2 at 850yum (for eight sources), and with IRAM/Gismo at 2 mm. The light and dark blue curves 
show the blackbody curves for optically thin and optically thick dust, respectively. Error bars are printed in black, and are smaller 
than the symbols showing the photometric data points. 


6. Signatures of gravitational lensing 

6.1. Flux densities and morphologies 

The first observational hint that our targets might indeed be 
strongly gravitationally lensed galaxies is their sheer brightness 
and small angular size below 18" (the FWHM of the SPIRE 
point spread function at the highest frequency, 250 /im), whereas 


90% of the 234 Planck high-z candidates with Herschel/SPIRE 
follow-up consist of multiple sources. Unlensed high-redshift 
galaxies with FIR flux densities as high as those of our sources 
are very rare. For example, for flux densities S 350 > 400 mJy at 
350/rni and near the pe ak of the FIR SEDs of most GEMS, the 
Bethermin et al. (2012) models predict 3 sources sr -1 at z>2 that 


are not, and 53 sources sr 1 that are gravitationally lensed. For 
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Fig. 5. Characterizing our sources by their FIR luminosities, 
dust temperatures, and redshifts. Top : FIR luminosities versus 
redshifts. The Dusty GEMS targets are shown as red stars on top 
of the comparison samples of strongly lensed sources taken from 
the literature (see |Greve et al.|2012| an d refe rences therein), lens 
candidates from the S PT survey ( Vieira et al.|2013 ), the Eyelash 
( [Ivison et~a l. 2010b), and the two Planck lenses that were dis¬ 
covered prior to our work fCombes et al.|2012} Fu et al.|2012| ). 
The solid line shows a fit to the general population, while the 
dotted lines show magnifications factors of 10 and 50. Bottom'. 
FIR luminosity as a function of redshift. In the lower panel, our 
sources fall near the most luminous gravitationally lensed galax¬ 
ies at similar redshifts; however, these are strongly dominated by 
AGN, like APM 08279+5255, whereas our sources are bona fide 
starbursts, at least in the image plane. This explains the greater 
luminosities at lower dust temperatures seen in the top panel. 


a given source at these flux densities, it is therefore much more 
likely to find a gravitationally lensed galaxy than a galaxy in the 
field. We emphasize that this is true only for galaxies with con¬ 
firmed high redshifts like the GEMS. The overall population of 
bright FIR galaxies is at much lower redshifts and has bluer FIR 
colors. 

Given that more typical high-z galaxies in the field have 
FIR flux densities of few lOsmJy, the observed flux densities 
of our sources make us expect magnification factors of at least 
a few and perhaps greater than 10, even if these were intrinsi¬ 
cally amongst the most luminous sources on the sky. For more 
moderate, more typical intrinsic flux densities, the gravitational 
magnification would be accordingly higher. The one alternative 
interpretation could be that these are multiple, highly concen¬ 


trated galaxies within projected distances of a few tens of kpc, so 
that they are not resolved into individual galaxies by the SPIRE 
beam, which corre spond s to about 120 kpc at z ~ 2 (for an exam¬ 
ple see Ivison et al. , 2013| ). 

More direct evidence that our sources indeed have a lensing 
nature comes from their morphologies. Strongly gravitationally 
lensed galaxies may either be single or multiple compact im¬ 
ages seen near the intervening foreground lensing source, giant 
arcs extending over several arcseconds, or even partial or com¬ 
plete Einstein rings. However, reaching the required spatial res¬ 
olutions of about 1-2" or better requires either deep optical or 
near-infrared imaging for sources that are not too heavily ob¬ 
scured, or, alternatively, submm or millimetre interferometry. 

Fig. [T] shows the Herschel image as contours overlaid on op¬ 
tical imaging, and Fig.[2]the dust morphology obtained with the 
SMA at 850yum. These figures clearly show that our sources are 
either single, compact objects at 2" resolution (as is the case 
for, e.g., PLCK_G244.8+54.9 and PLCK_G138.6+62.0), or gi¬ 
ant arcs (e.g., PLCK_G145.2+50.9 and PLCK_G080.2+49.8). 
PLCK_G244.8+54.9 and PLCK_G138.6+62.0 have flux densi¬ 
ties of 1054+10 and 664+8 mJy, respectively, making it implau¬ 
sible that these are individual high-redshift galaxies if not bene- 
fitting from a strong boost from gravitational lensing through a 
massive foreground source. At a 2" beam, intrinsic source sizes 
of unlensed galaxies would be at most 16 kpc, making it im¬ 
plausible to see a small group of high-redshift galaxies within a 
single beam. 

Our sources are associated with overdensities of massive in¬ 
tervening galaxies at intermediate redshift, either galaxy groups 
or clusters. Table |T] lists the redshifts of the brightest intervening 
galaxies along the line of sight to our targets, which were taken 
from the SDSS. They are either photometric redshifts or (in a 
few individual cases highlighted in Table [I]), spectroscopic red¬ 
shifts. A detailed analysis of the lensing structure based on our 
own proprietary optical/NIR photometry obtained at the CFHT 
and the VLT is on-going. The only target associated with a fore¬ 
ground object for which we currently have no good redshift es¬ 
timate is G145.2+50.9, which forms a near-complete Einstein 
ring with a diameter of 10" around a very red foreground source, 
which we currently only detect in the J-band. For this galaxy, our 
optical/NIR photometry is not yet complete, and we must there¬ 
fore defer a detailed analysis of the foreground source to a later 
publication. 

In the rest of this section we will present additional empirical 
evidence to illustrate that all of our sample consists of strongly 
gravitationally lensed galaxies, including those for which we do 
not yet have morphological constraints. 

6.2. FIR dust continuum 

Greve et al. ( |2012| ) suggested the use of the relationship be¬ 
tween and Lfir as an empirical indicator of whether a source 
is strongly gravitationally lensed. According to the Stefan- 
Boltzmann law, the luminosity emitted by a modified black- 
body at a given temperature depends only on the size of the 
emitting surface and the emissivity. Therefore, if our galaxies 
form a subset of the generic population of dusty, intensely star¬ 
forming high-redshift galaxies, the gravitational magnification 
factor should lead to an apparent increase in dust luminosity at a 
given temperature for gravitationally lensed galaxies, compared 
to unlensed galaxies. 

In Fig. [5] we show where our sources fall on a plot of dust 
temperature versus FIR luminosity, compared to the galaxies of 
Greve et al] ( |2Q12| ) and other sources from the literature. It is 
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immediately clear that our sources are significantly brighter at 
a given temperature than galaxies in the field, and even brighter 
than most of the previously known gravitationally lensed galax¬ 
ies. The dashed lines illustrate luminosities that are 10 times and 
50 times brighter than expected from a simple least square fit 
to the population of field g alaxies (green solid lin e). We high¬ 
light the Cosmic Eyelash of Swinb ank et al.| ( |2010| ) and the two 
previously confirmed gravitationally lensed galaxies with Planck 
counterparts et al. 2012[|Combes et al.|2012| ). 

For a more explicit, although not necessarily very precise es¬ 
timate that exploits this relationship between dust temperature, 
luminosity, and source size, we can use the Stefan-Boltzmann 
law directly to infer the magnification factor from the observed 
luminosity and dust te mpera ture. The approach is described in 
detail in Greve et al. ( |2012| . Given the somewhat higher dust 
temperatures of our lens candidates compared to the SPT sources 
of Greve et al. ( 2012 ), we have a larger range of “effective” 
Stefan-Boltzmann constants (which correct for the lower emis- 
sivity of a modified relative to a genuine blackbody), <x e ff/<x = 
0.5-0.65, and apparent emitting surfaces that are between 8 and 
40 times larger than those of typical FIR galaxies, implying mag¬ 
nification factors that are roughly similar. These results are also 
shown in Fig. [5] We note that the estimates are very uncertain, 
since they assume that each high-redshift galaxy has a similar 
emissivity, that the dust becomes optically thick at similar wave¬ 
lengths, and that the interstellar medium (ISM) of the galaxies 
is dominated by a uniform dust component with a single tem¬ 
perature. Each of these assumptions can only be approximately 
true, and this will lead to uncertainties of at least factors of a few 
in this estimate. This is illustrated by the magnification factor of 
about 10 of the “Cosmic Eyelash,” which we constrained with 
the same approach, and which falls into a very similar region of 
the diagram as our sources. The “Cosmic Eyelash” has intrinsic 
magnification factors of 20-60 ( [Swinbank et al.[20 11), with an 
average factor of 32.5. In spite of these caveats, Fig. [5] provides 
additional indirect support that our sources are indeed strongly 
gravitationally lensed galaxies. 


6.3. Molecular gas lines 

More empirical evidence for strong gravitational magnification 
of our galaxies comes from the emission-line properties that we 
measured with EMIR on the IRAM 30-m telescope. To demon¬ 
strate that their galaxies discovered through the H-ATLAS sur¬ 
vey are indeed strongly gravitationally lensed, [Harris et al. 
(2012) compiled the CO(l-O) line luminosities U and FWHMs 
of high-redshift galaxies, finding that typical galaxies follow a 
broad trend between luminosity and line width (Fig. 6]). Their 
own galaxies and other strongly gravitationally lensed galaxies 
from the literature stand out by 1-2 orders of magnitude above 
this relationship, owing to the boost in line luminosity by the 
gravitational lens. In contrast, lines are expected to be as broad 
as in field galaxies (or narrower, because smaller regions of the 
large-scale velocity field of the galaxy are being observed). If 
a small region in a galaxy is magnified by strong gravitational 
lensing, then only parts of the rotation curve will be sampled, and 
turbulent motion will also be smaller. Therefore, observed line 
widths should be narrower than if galaxy-wide radii are probed, 
as in field galaxies. 

Fig. [6} inspired by Figure 7 of Harris et al. ( |2012| ), shows 
where our sources fall relative to the Harris et al. sample and 
those in the literature. Our sources span a large range in FWHM 
in this diagram, but are all clearly within the regime of gravita¬ 
tionally lensed galaxies. The FWHMs of some of our targets may 
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Fig. 6. Top: Order-of-magnitude estimate of the gravitational 
magnification factor derived from the dust temperature and FIR 
luminosity of our sources compared to FIR/submm galaxies in 
the field (see Sect. |6.2| for details). Bottom : CO emission-line 
luminosity as a function of FWHM line width, using the same 
symbols as in the top panel, with additional l ight blue dots show¬ 
ing those sources of |Bothwell et al.| ( |2013| ), which are neither 
gravitationally lensed nor AGN hosts. The solid and dashed lines 
show the best linear fit to the sample of field galaxies, as well as 
shifts of this line corresponding to magnification factors of 5, 
10, and 20. Error bars show the ler uncertainties of our single¬ 
component Gaussian line fits. 


be as large as those of overdensities of multiple lensed sources 
that are at best_moderately gravitationally lensed like the exam¬ 
ple in Ivison et al. |2013[ ). Their luminosities are factors of a few 
brighter, placing them firmly above the field galaxies, and into 
the lens regime. The example of the Cosmic Eyelash, which has 
a magnification factor of 20-60 shows however that the lens¬ 
ing factor derived from the position of a source in this diagram 
can only be an approximate indicator. The line luminosities and 
profiles provide additional evidence in support for the lensing 
hypothesis. 


7. Dust heating and AGN content 

Many du sty high-z starburst ga laxies host powerful AGN. For 
example Alexander et al. ([2005) find that 75 % of submillimetre 
galaxies have AGN detected at X-ray wavelengths (but see also 
Laird et a l.|2010 who find significantly lower fractions of 10- 


20 %). Powerful, obscured AGN could contribute significantly 
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to Lfir, which would lower the magnification factors that are 
needed to explain the extraordinary apparent brightness of our 
targets. The inferred dust temperatures of our sources are akin 
to those of other high-z starburst galaxies and would be very 
low compared to those of obscured quasars, which are typically 
above 70 K, often far above. However, our estimates may be bi¬ 
ased towards low temperatures because of the absence of mid- 
infrared constraints shortward of 50-80yum in our data sets. This 
makes it worthwhile to investigate in more depth what the possi¬ 
ble AGN contribution might be to the infrared luminosity budget 
of our targets. 

Not all AGN play a significant role for the total infrared lu¬ 
minosity budget of their host galaxies. Ideally, searching for ob¬ 
scured AGN activity in our targets would require PACS or at 
least deep Spitzer 24-pm photometry, to sample the dust SED 
shortward of about 100 /am in the rest-frame. Since our sam¬ 
ple was only discovered with SPIRE in the last weeks of the 
Herschel mission, such data sets are not available. Additionally, 
in spite of their exceptional brightness (for high-redshift targets), 
they are still too faint to obtain robust mid-IR photometric con¬ 
straints with SOFIA. We might therefore not be able to give 
tight constraints on the presence of weak AGN activity in all 
of our sources, however, most important for our purposes here 
is to show that AGN radiation is not energetically dominant 
in our sources (i.e., they do not contribute > 0.5 L^/), which 
would have important implications for the estimated star forma¬ 
tion rates, plausible ranges of intrinsic infrared luminosities, and 
gas excitation. 


7.1. WISE 22pm photometry 

The Wide-Field Infrared Survey Explorer (WISE) in the 22-pm 
band covers rest-frame wavelengths between 4 and 7 pm for our 
sources, i.e., the far blue tail of hot dust emission from powerful 
AGN. Consequently, WISE has been successfully used to iden¬ 
tify heavily obscured AGN at high-z (e.g., |Yan et al.|2013) . All 
of our sources are bright enough to be detected in WISE imaging 
at 22pm, and seven are included in the WISE catalogues. The 
remaining four can be seen as very faint sources in the WISE 
images, but are not bright enough to obtain robust flux measure¬ 
ments. We will in the following use the flux measurements from 
the WISE catalogs and treat the 5<x rms of this catalog at 5.4 mJy 
as upper limit for the fainter sources. Individual flux densities are 
listed for all sources in Table Q] 

Although we do detect nearby bluer sources in each case at 
3 A pm, 4.6 pm, and 12 pm, the spectral energy distribution in 
all but one source falls off too steeply in the three blue WISE 
bands to contribute significantly to the 22pm detections. Their 
colours suggest that these are low-to-intermediate-redshift ob¬ 
jects, which we consider part of the intervening lensing struc¬ 
ture. 


7.2. Constraining AGN dust heating with DecompIR 

We used the publicly available package DecompIR ( [Mullaney] 
et al.1 [2011| ) to constrain the potential contribution from AGN 


emission to the infrared SED, using our six available data points 
between 22pm and 2 mm, as listed in Table [T] DecompIR per¬ 
forms x 1 fits to infrared SEDs using composites of empirically 
constructed templates of starburst galaxies and AGN. We only 
have six photometric data points, but we sample the Rayleigh- 
Jeans tail and the peak of the spectral energy distribution of the 
emission from the coldest dust component, and therefore con¬ 


strain the most abundant dust component well. We have spec¬ 
troscopic redshifts and hence can estimate robust dust temper¬ 
atures, and therefore we consider this part of the SED to be 
well constrained. Our main uncertainties are blueward of the 
dust peak, which we only sample with the 22 pm data point from 
WISE and the 250 pm data point from SPIRE. The 250 pm band 
corresponds to 50-80 pm in the rest-frame for our targets, is 
slightly blueward of the peak of the dust emission in most of our 
sources (Fig. [4]), and falls near the expected peak of the infrared 
emission from hot dust in the AGN torus in most templates (e.g., 
|Polletta et al.|2007(|Nenkova et al.|2008] ). The SEDs of power¬ 
ful AGN host galaxies at high redshifts are t ypically dominated 
by AGN emission a t these wavelengths (e.g., Saji na et al.|2012| 
prouart et al.|2014| ). 

The DecompIR software package provides a small set of SED 
templates, derived from nearby starburst galaxies. One compli¬ 
cation in using these templates for our sources is that we find a 
range of dust temperatures, which are not all well represented by 
existing templates. This would lead to considerable (but some¬ 
what artificial) discrepancies in the template fitting that would 
be hard to overcome without a self-consistent model of the mid- 
to-far-infrared SED of high-redshift galaxies, which does not yet 
exist. To avoid such mismatches, while make optimal use of our 
existing constraints without overinterpreting the mid-IR data, we 
therefore constructed a simple starburst template for each of our 
galaxies from the “SB2” template of DecompIR (corresponding 
to the spectral energy distribution of NGC7252, |Mullaney et aL] 
2011| ), but correcting for the mismatch in dust temperatures. 


We therefore fitted and removed from the template the modified 
blackbody contribution from cold dust in the FIR. The residual 
is a template of the SED in the mid-infrared only, to which we 
added the modified blackbody emission obtained from our FIR 
fits, using the measured temperature for each source. We selected 
this particular template, because it provided the best match to 
the FIR and sub-millimeter measurements of the dust peak and 
Rayleigh-Jeans tail, with no regard of the goodness of fit of the 
22 pm observations. For the AGN component, we simply used 
the |Polletta et al.|p007| ) type-2 QSO template which is already 
part of DecompIR. 

Fig. [7] shows that we can rule out a bolometrically dominant 
AGN contribution from buried quasars, in all cases. Only three 
galaxies have best-fit results formally inconsistent with pure star- 
bursts. In PLCK_G092.5+42.9 and PLCK_G244.8+54.9, the pu¬ 
tative AGN contribution to the FIR luminosity is below 10%, 
and it is below 30% in PLCK_G145.2+50.9. In all other cases, 
DecompIR finds the best fit for a pure starburst SED. 

Fig. [7] also shows the 90% completeness limits of the IRAS 
faint source survey (Moshir et al. |1992[ ), 120 mJy and 440 mJy 
at 60 pm and 100 pm, respectively (green downward arrows). 
Although these limits are not constraining for the SED fits that 
we obtained with SPIRE, SCUBA-2, GISMO, and WISE (and 
were therefore not included in our DecompIR fits) they do illus¬ 
trate that a hypothetical bright quasar component would have led 
to a 60 pm detection at least in our brightest targets. In Fig. [8] we 
illustrate the impact of fiducial AGN contribution of 0.1, 0.3, and 
0.5 xLfir on the example of G138.6+62.0, a source with rather 
average FIR brightness and 250-to-350 pm color in our sample. 

Obviously, these fits are uncertain, given that they cover the 
far blue wing of the dust SED expected from bright AGN, and 
a spectral range that has important contributions from rich mid- 
infrared spectral features, in particular from Polycyclic Aromatic 
Hydrocarbons (PAH), which dominate the infrared spectral en¬ 
ergy distributi ons o f starburst galaxies between about 3pm and 
10/rm (e.g., [Armus et al.|2007| ). 
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Fig. 7. Infrared-to-millimetre spectral energy distributions of all of our 11 sources obtained using DecompIR (Mullaney et al. 


2011 ) 


Red dots show our data points, where the error bars along the abscissa indicate the width of each band. The error bars along the 
ordinate are often smaller than the symbol size. Green downward arrows show the 90% completeness limit of the IRAS all-sky 
survey at 60 and 100 ji m, respectively. Red downward arrows at 22/dm show the 5 cr upper limits provided by the WISE catalogs, 
where the counterparts were fainter than the 5.4 mJy flux limit of the WISE catalog. Black, blue, and yellow lines show the starburst 
and AGN component (if an AGN was fitted), and the sum of both, respectively. 


The equivalent widths of the PAH bands have a wide 
dispersion in high-redshift galaxies (e.g., Fiolet et al.| |2010| ), 
and are generally larger than those measured in low-redshift 
ULIRGs, including NGC 7252, which we used as template in 
our DecompIR fits. Additional uncertainties come from the stel¬ 
lar continuum, which reaches comparable strength to the dust 
continuum in the wavelength range covered by the 22 ji m ob¬ 
servations, about 5-7 /tm, and whose shape depends strongly 
on the specific star-formation rate, star-formation history, and 
geometry of dust obscuration in the host galaxy. Compared to 
the average mid-infrared SED of starburst galaxies ( [Brandi et al.| 
200 6]), the stellar contin uum in NGC 7252 seems relatively faint, 
and Sajina et ah] ([20 12 ) also found that the mid-IR continuum of 
high-redshift starburst galaxies seems to be brighter than in low- 
redshift galaxies. Other spectral features might also contribute 


to boosting the observed 22 ji m flux density. For example [Fiolet 
|et al.| ( [20 10| ) report the detection of warm H 2 line emission at 
6.9 /a m in the rest-frame. Moreover, the foreground galaxies may 
also contribute to the 22 /dm flux, given the 12" beam of WISE. 
However, all these uncertainties would only act to lower the pu¬ 


tative intrinsic contribution of AGN heating relative to what we 
observe. It will of course be very interesting to study the mid- 
infrared spectral properties of our sources in depth once that 
JWST/MIRI will become available. 

Relying upon these arguments, we conclude that the dust 
heating in our sources is not dominated by heavily obscured 
quasars, but that the Planck Dusty GEMS are genuine starburst 
galaxies with at most a minor contribution of AGN heating to the 
dust emission and the bolometric energy budget of our sources. 
Of course this only refers to the regions that we are seeing under 
the gravitational microscope of the foreground lensing systems. 
The global energy budgets of AGN and star formation in these 
galaxies may be different, if the AGN lies alon g a line of sight 
that is not strongly magnified (see also |Serjeant|2012[ ). 

8. FIR-radio correlation 

Given the extraordinary brightness of our targets on the sky, we 
can use the 1.4 GHz VLA survey of the northern sky, FIRST 
( Becker et al.|1995| , to search for counterparts of our sources at 
an observed frequency of 1.4 GHz (between 4.5 and 6.4 GHz in 
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Fig. 8. A schematic illustrating the the impact of 0.1, 0.3, and 
0.5 x Lfir of AGN contamination to the shape of the SED of 
our sources, on the example of G138.6+62.0 and G 165.7+67.0. 


the rest-frame). We find that six of our sources have counterparts 
within 5" of a source in the FIRST catalogue, and the other five 
have sensitive upper limits (see Table [4]). Here 5" corresponds 
to the beam size of FIRST. We assume a radio spectral index 
a i .4 = -0.8 (as is appropriate for star-forming regions, |Condon| 
\1992} to convert these flux densities to a monochromatic radio 
flux density at rest-frame 1.4 GHz, as well as to estimate a ra¬ 
dio power. For galaxies without detections we use the 3 cr upper 
limits implied by the rms given in the FIRST catalogue. Results 
for individual sources are listed in Table 0] All sources are listed 
in the FIRST catalog as marginally resolved, with deconvolved 
sizes between 1.5" and 7". The 5" beam of FIRST is not al¬ 
ways sufficiently small to rule out contamination from a radio 
nucleus in the foreground lensing galaxy. However, most sources 
are isolated enough to conclude that the radio counterpart is at 
the position of the high-redshift galaxy, and not coincident with 
a massive, intermediate-redshift galaxy that could host a radio- 
loud AGN. This is particularly the case for PFCK_G244.8+54.9. 

The far-infrared radio correlation is commonly parametrized 
as q = logLi R /(3.75 X 10 12 W) - logLi. 4 /(WHz _1 ), where L lR 
is measured in the range 8 - 1000 /tm in the rest-frame, and the 
centimetre continuum at a rest-frame frequency of 1.4 GHz. 
Early He rsche l results suggest ( q ) = 2.4 + 0.12 in the local 
Universe ( [Jarvis et al. 2010), although values can change con 


siderably even within individual galaxies (e.g., Tabatabaei et al. 
201 3). Towards hi gher redshifts, results show less of a consen¬ 
sus. |Kovacs et al.| (2006|), |Vlahakis et~aTT1 ( |2007| ), |MichaIowski| 


et al.| ( |2010| ), and [Bourne et al. ( 201 1) report an increase in ra¬ 


dio power at a given FIR luminosity at z ^ 1, with ( q ) = 2.0. 
However, several studies (Ivis on et al. 2010a|c[ Sargent et al. 


Table 4. Radio continuum properties extracted from the FIRST 
catalogue. The columns are: source name; relative distance be¬ 
tween the SPIRE position at 250/rm and the position of the 
FIRST counterpart, if detected within 5"; integrated flux density 
at 1.4 GHz from the FIRST catalogue and gravitational magni¬ 
fication factor of the radio component, ji 1 . 4 ; radio luminosity at 
1.4 GHz in the rest-frame and gravitational magnification factor 
of that component, /ti. 4 , for sources without 1.4 GHz detections, 
the luminosities are derived from 3x the rms given in column 3; 
the Ratio between radio and FIR luminosity (see Sect. [5]), i.e. the 
^-parameter. 


Source 

Dist. 

[arcsec] 

MlA S\ n 4 
[mJy] 

M 1.4 7-1.4 
[10 25 WHz -1 

q 

] 

PLCK_G045.1+61.1 . . . 


<0.135 

<3.2 

>2.4 

PLCK_G080.2+49.8 . . . 


<0.135 

<1.8 

>2.4 

PLCK_G092.5+42.9 . . . 

0.9 

1.50+0.16 

11.2 

2.35+0.11 

PLCK_G102.1+53.6. . . 


<0.149 

<2.5 

>2.5 

PLCK_G113.7+61.0. . . 

3.0 

1.9+0.14 

6.9 

2.2+0.08 

PLCK_G138.6+62.0 . . . 

2.7 

2.01+0.16 

7.3 

2.1+0.08 

PLCK_G145.2+50.9. . . 


<0.144 

<3.9 

>2.8 

PLCK_G165.7+67.0 . . . 

4.1 

3.41+0.15 

10.1 

2.0+0.04 

PLCK_G200.6+46.1 . . . 

1.1 

1.23+0.14 

6.9 

1.9+0.11 

PLCK.G231.3+72.2. . . 


<0.151 

<2.5 

>2.5 

PLCK_G244.8+54.9. . . 

3.2 

2.26+0.14 

13.6 

2.3+0.06 


2010; Thoms on et al. 201 4[ e.g.,) find average values in the range 


q = 2.4-2.7, out to z ~ 2. Most of these studies use stacking anal¬ 
yses to overcome the inherent faintness of high-redshift galaxies 
in the FIR/submm and in the radio, and are based on ground- 
based data that do not sample the FIR dust SED very well. In the 
Cosmic Eyelash, a single, strongly lensed galaxy at z = 2.4 with 
Herschel-SPIRE coverage, Ivison et al. ( |2010b| find q = 2.4. 

In Fig. [9] we show where our sources fall relative to the local 
FIR-radio correlation, assuming that differential lensing between 
the FIR and radio synchrotron emission does not play a major 
role in our case. Four of our sources are brighter by up to about 
0.4 dex in the radio, and fall outside the lcr scatter around the 
distribution of q parameters of low-redshift galaxies. Six galax¬ 
ies fall within lcr from this relationship, however, three of these 
galaxies only have upper limits on q , and would fall outside the 
1 cr scatter if we had adopted lcr, instead of 3 cr upper limits. One 
galaxy has q>2.7 at 3 cr. 

If the 1.4 GHz emission is caused by star formation, then the 
lower q values cannot be caused by our assumption of the spec¬ 
tral index or 1.4 = -0.8. Only for flat radio spectral energy dis¬ 
tributions (aq .4 = - 0.2 or above, which would only be possible 
for an AGN) would our highest-^ sources be consistent with the 
FIR-correlation and q = 2.4. If this were the case, our sources 
with the highest q parameters (i.e., the brightest sources in the 
FIR relative to the radio) would fall well above q = 3. 

One way of obtaining lower g-values , as has been pointed 
out several times before (e.g., Vlahakis et al.|2007} Bourne et al. 


2011) is through radio emission from a central AGN. This has 


an interesting theoretical aspect, since it has been proposed (e.g., 
Silk & Norman 2009) that gas compression caused by expand¬ 


ing radio sources in massive high-z galaxies may in fact con¬ 
tribute to boosting the star-formation efficiency in intense high- 
z starbursts. However, although contamination with nuclear ra¬ 
dio emission may be a particular concern for stacking analyses 
(where one or a few relatively bright radio sources may affect 
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|Lehnert et al.|20Q9t |Swinbank et al.|2011| ), which enhances the 

scale height of the gas, and lowers their volume density. In this 
case, energy loss of cosmic rays through synchrotron radiation, 
bremsstrahlung, and other processes could be either enhanced or 
diminished, depending on the local magnetic fields and density 


Fig. 9. Far-infrared radio correlation of our sources, as discussed 
in Sect. [ 8 ] /jfir and /i 1.4 correspond to the magnification fac¬ 
tors appropriate for the dusty and 1.4 GHz synchrotron emis¬ 
sion, respectively. Solid blue and red lines show the relationship 
for q = 2.4 and 2.0, respectively, while the red dotted line if for 
q = 2.7. The dot-dashed lines_show the ±1 cr scatter about q = 2.4 
at low redshift from Jarvi s et a l.]( [2010| . 


the average result of the overall sample), in our case this expla¬ 
nation seems to be less likely. First, the radio sources are spa¬ 
tially resolved in FIRST, suggesting they are extended, just like 
the star formation. Even if the Dusty GEMS have extended radio 
sources, the gravitational lensing would need to amplify the ra¬ 
dio emission by similar amounts to that of the star formation in 
most of our sources; this appears relatively unlikely, given that 
the AGN radio morphologies of high-z galaxies are very differ¬ 
ent from those of the dust and stellar components (e.g. |Sajina| 
let al.|2007| ). 

As a second counter argument, our SED fits have already dis¬ 
favoured the presence of very luminous AGN. Although galax¬ 
ies are known to host central radio sources without bright bolo- 
metric emission, the host galaxies of such AGN tend to have 
little dust, gas, and on-going star formation. And although rela¬ 
tively faint AGN may be present without leaving strong signa¬ 
tures in the dust SEDs, only a subset of those should be radio 
loud. Even optimistic results (e.g., from Sajina et al.p 007[ ), im¬ 
ply that only up to about 30 % of dusty high-z starburst galaxies 
may host moderately bright (Pi .4 « 10 25 WHz -1 , comparable to 
our sources for magnification factors of order 10 ) nuclear radio 
sources that are as bright or brighter than the star formation it¬ 
self. In addition, these sources are classified as AGN in the mid- 
infrared. We note that we have not found any obvious trends 
between the ^-parameter and possible probes of AGN contami¬ 
nation like dust temperature or Lfir- We also do not find a cor¬ 
relation between the flux density at 22 /am and the 1.4 GHz radio 
emission. 

As an alternative, the offsets towards larger and smaller q- 
parameters may be caused by the star -forming environments 
themselves. |Lacki & Thompson] p010| ) suggest that enhanced 
synchrotron emission from cosmic rays in star-forming galaxies 
at high redshift could be one outcome of the strong turbulence 
observed in these galaxies (e.g., Forster Schreiber et al.| |2009; 


distribution of the ISM; this could either decrease or increase 
the energy losses of the cosmic rays, making our sources either 
brighter or fainter (e.g., [Murphy[2009] ). Other possible explana¬ 
tions, which might also involve the sources with unusually high 
^-parameters, include evolution in the dust properties, the age 
of the starburst (the radio should not probe starbursts with ages 
of less than a few times 10 7 years or greater than a few times 
10 8 years, because of the timing of core-collapse supemovae, 
Br essan e t al. 2002| ), and a top-heavy initial mass function (e.g., 
Baugh et al.|2005| ). It will be interesting to obtain high-resolution 
centimetre-wave maps of these sources, along with our gas and 
dust interferometry, to constrain the potential impact of differen¬ 
tial lensing and to further elucidate how the ^-parameter depends 
locally on the ISM properties in our sources. 

9. Gas masses and gas-to-dust ratios 

To determine the molecular gas mass (which is dominated by 
H 2 ) from the CO luminosity, we have to assume an empiri¬ 
cal conversion factor, which is notoriously difficult to justify 
from first principles, and which is therefore still heavily de¬ 
bated inAhe_^iterature, in particular for high-redshift galaxies 
(e. g.,|Daddi et al.|20T0||Genzel et al.| 2010 | Glover & Mac Low 

2011; Naraya nan et al.||2012| ). The canonical value adopted for 

most high-z galaxie s is a & O.SMQ/^Kkms " 1 p c 2 ). This value 
was first derived by |Downes & Solomon] ( |1998| ) for the dense, 
circumnuclear, optically thick molecular gas discs in nearby 
ULIRGs, and is commonly adopted also for dusty starburst 
galaxies in the early Universe. 


However, several studies in recent years (starting with Daddi 
et al. |2010| and |Genzel et al.pOlO] ) have called into question 
whether a single CO(l- 0 )-to-H 2 conversion factor may be ap¬ 
propriate to use for all high-z galaxies. Several attempts have 
therefore been undertaken to constrain ac o either from dynami¬ 
cal mass estimates, or on theoretical grounds. One empirical ap¬ 
proach suited for galaxies with well constrained dust mass mea¬ 
surements like ours, is to assume that high-z galaxies fall onto a 
similar relationship between gas-to-dust mass ratios and metal- 
licities as found for SINGs galax ies in the local Universe ( [Leroy 
et al.|2008] ). |Magdis et al.| ( [2011| ) used this approach to confirm 
their aco determinations, which they previously obtained from 

dynamical mass estimates. _ 

We can use our measurements of the dust mass, ( Sect . |6.2| 
and Table [2]), and the CO line luminosity, L' (Sect. |4.3| and 
Table [3), to esti mate ratios of L' /Md, which scale wi th gas- 
phase m etallicity ([Leroy et al.|2 008; Ma gdis et al.|2011| ). Using 
Lig. 3 of Ma gdis et al.| ( |2011| ), we find that our measured range 
L'/Mj = 40- 140 corresponds to gas-phase metallicities 12 + 
log(Q/H ) ~ 8.9 - 9.3. These values are appropriate for J3 = 2.0 


(Sect. 6.2). Lurthermore, we adopted a ratio r 32 /io = 1 to convert 
from the luminosities of the observed mid-J CO lines to CO(l- 
0). Thi s factor is expected for optically thick gas (e.g., 

& Vanden Bo ut|2005| . In the Cosmic Eyelash, Danie 
j2011) measured r 32 /io = 0.7. 


Solomon 


son et al. 


High gas-phase metallicities correspond to small values of 
aco of < 1.0. Lor example, if we use the linear fit between 
aco and metallicity of Genzel et ah] ( |2012| ), we find conver¬ 
sion factors of about 0.4 at face value, although with large 
uncertainties. Likewise, Lig. 3 of [Magdis et al.| ( |2011| sug¬ 
gests aco < 1.0M o /(Kkms _1 pc 2 ) for galaxies with L'/Afj 
ratios and metallicities as high as in our sources. This sug¬ 
gests that using the ULIRG conversion factor of aco = 
0 . 8 M o /(Kkms _1 pc 2 ) is more appropriate than much higher 
factors of 3 - 5 M 0 /(Kkms _1 pc 2 ), as previously adopted for 
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more moderately star forming, disk-like high-redshift galaxies 
and the Milky Way. We stress that these results are measured in 
small regions of high-z galaxies, and are not necessarily repre¬ 
sentative of the average values in these galaxies ( Serjeant|2012 ). 

When using a common conversion factor aco ~ 
0.8M o /(Kkms _1 pc 2 ), we find molecular gas masses of 
2-lOxlO 1 V _1 M 0 for Planck's Dusty GEMS (including the 
gravitational magnification factor /d). Results for individual 
galaxies are listed in Table [3] 


10. Integrated star-formation law 

The tight correlation between molecular gas mass surface den¬ 
sity and star-formation intensity over scales of around 100 pc 
to entire galaxies highlights how star formation depends on the 
avail able molecular gas reservoirs ( i.e., the Schmidt-Kennicutt 
law, Schmidt 1959; Kennicutt 1998 ), although there is no con¬ 
sensus about a unique physical mechanism putting this relation¬ 
ship in place. 

However, starting with, e.g., |Daddi et al. ( 2010| ) and Genzel 


[eTai: (2dT0|), several studies have emphasized in recent years 
that not all high-redshift galaxies may strictly obey the same 
empirical star-formation law, but that at a given gas surface den¬ 
sity, starburst galaxies may be more efficient in turning their gas 
into stars. The reasons for this are still unclear, with possibilities 
ranging from changes in the stellar initial mass function ( |Baugh| 
|et al.|p005] ) to different star-formation efficiencies. For exam¬ 
ple, there may be more star-forming clouds in the ISM than in 
more quiescent high-z galaxies ( [Lehnert et al.|201 3) , or the star- 
formation efficiency per free-fall time could be higher ([pekel 


et al.|2009| . What we can confidently assert (given the absence 


of dust with T»50K, faint 22-yum flux density, or excess radio 
emission) is that AGN do not play a dominant role in boost¬ 
ing the FIR continuum in the strongly amplified regions that we 
are seeing in our galaxies (Sect. [7]). Given the close astrophysi- 
cal connection between dust and gas in star-forming regions, we 
see no reason to believe that differential lensing plays an impor¬ 
tant role for our analysis of the gas conditions in the regions, 
which are amplified by the gravitational lenses. Of course, we 
must keep in mind that the gas and dust measurements obtained 
in small, selectively amplified regions of these galaxies are not 
necessarily representative of the average gas an d dust properties 
in these galaxies on global scales. ( |Serjeant|2012| ). 

The relationship of far-infrared luminosity with CO line lu¬ 
minosity is a simple, empirical, and robust way of investigating 
this “star-formation law.” It is directly related to the integral form 
of the Schmidt-Kennicutt law, because the far-infrared luminos¬ 
ity probing the emission from warm dust is an excellent tracer 
of star formation. CO line emission is mainly excited through 
collisions with H 2 and is the most commonly adopted proxy of 
the total molecular gas mass (e.g., |Qmont|2007| ). 

In order to investigate whether our sources are representa¬ 
tive of star formation in high-z starburst galaxies or more ordi¬ 
nary, but nonetheless very intensely star-forming galaxies on the 
“main sequence” (e.g., |Elbaz et aT7||201 1| ), we use the molecu¬ 
lar gas masses estimated from the CO luminosities (Mm) and 
far-infrared luminosities (Lrir), shown in Fi g. [TO] and com pare 
with the correlations found by Da ddi et al.| ( [20T0| ) and |GenzeI| 
|et al.| ( [20T0| ) for starburst galaxies, together with the more grad¬ 
ually, but nonetheless intensely star-forming high-z galaxies on 
the “main sequence.” For this comparison, we use the inte¬ 
grated star-formation rates and molecular gas masses, adopting 
aco = 0.8M o /(Kkms _1 pc 2 ) as justified in Sect. [ 9 ] 




9.5 10.0 10.5 11.0 11.5 12.0 

log IM [M ] 

3 H2 L surf 



Fig. 10. Left: Infrared luminosities (Lir) as a function of molec¬ 
ular gas mass for our sources. Red, yellow, and purple stars in¬ 
dicate the measured values (uncorrected for gravitational am¬ 
plification), and the same values corrected by fiducial lensing 
factors of 10 and 50, respectively. In agreement with the dust-to¬ 
gas ratios discussed in Sect. [9] we adopt the classical “UFIRG” 
conversion factor aco = 0.8M o /(Kkms _1 pc 2 ) for all sources. 
The blue and red line indicate the relationships forstarbursts and 
“main-sequence” galaxies emphasized by [Daddi et al.| ((2010 ) 
Small dots show the low- and high-z samples of |Daddi et al 


( 2010| ) for comparison. Fight blue dots along the upper line show 
low- and high-z UFIRGs and SMGs, respectively. The filled dark 
green dots (with error bars) are BzK galaxies at z ~ 2; for details 
of these samples see Daddi et ak| ( |2010| ). All error bars include 
the measurement uncertainties. The black cross in the lower right 
comer of the plot shows a fiducial factor 2 in statistical uncer¬ 
tainties. 


We find that all sources fall closer to the the upper “starburst” 
than the lower “main sequence” lines in Fig. [TO] To illustrate the 
effects of lensing, i.e., how these sources might appear in the im¬ 
age plane, we plot not only the observed luminosities (red filled 
stars), but also the luminosities for fiducial magnification factors 
Id = 10 and id = 50 (orange and light blue empty stars in Fig. [TO]). 
Even without the detailed lens modelling their position in the 
diagram relative to typical starburst or main-sequence galaxies 
does not depend sensitively on the precise magnification factor. 

We stress that in the present study, we only show inte¬ 
grated measurements, whereas the relationship that underlies the 
Schmidt-Kennicutt diagram is between the surface density of 
molecular gas mass and the star-formation rate. Our on-going 
interferometric follow-up programme will enable us to derive 
more detailed constraints, and to investigate whether the “star- 
burst” and “main sequence” classifications are unique for each 
source, or whether each individual Dusty GEM will show a 
range in star-formation efficiency. 

FIR luminosities and gas mass estimates also provide sim¬ 
ple, rough constraints on the gas depletion timescale, t& v \ = 
M g a S /SFR, and hence the time during which the current star 
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formation intensity can be maintained without replenishing the 
molecular gas reservoirs. Assuming that star-formation rates are 
constant, we find short gas depletion timescales, t & p i = (0.5-6) x 
10 7 years, significantly less than the typical stellar age of a far-IR 
or submm galaxy (a few times 10 8 years, Smail et al. 2004[ Lapi 


|et al.|[20l2| ), and also somewhat, but not dramatically, shorter 
than the gas depletion timescal es found for unlensed submm 
galaxies (e.g., |Greve et al.|2005[ ). This highlights again the fact 
that our sources have all the hallmarks of being “ordinary” dusty 
starburst galaxies placed under particularly powerful cosmic mi¬ 
croscopes. 


11, Summary 

We have presented a first analysis of an extensive multi¬ 
wavelength follow-up campaign of a new sample of the bright¬ 
est high-redshift FIR/submm galaxies, discovered through the 
unique synergy of the Planck and Herschel satellites. Planck's 
all-sky nature and multi-frequency coverage allows us to se¬ 
lect rare peaks in the submm background and Herschel ob¬ 
servations lead to sub-sample of strongly lensed candidates - 
Planck's Dusty GEMS. Their FIR peak flux densities are up to 
S' 350 = 1130 mJy at 350/tm, including six sources that are above 
the completeness limit of Planck at the highest frequencies. Our 
sample extends the very successful searches for gravitationally 
lensed high-z galaxies already carried out with Herschel and the 
SPT towards the brightest, rarest targets on the FIR/submm sky, 
which emphasized the need for a genuine all-sky survey to sys¬ 
tematically probe such exceedingly uncommon sources. 

All sources in our sample are bright, isolated point sources in 
SPIRE 250-yum maps (18" FWHM), and have the typical FIR- 
to-mm SEDs of dusty, intensely star-forming galaxies at high 
redshift. They have redshifts in the range z = 2.2-3.6, based on 
multiple bright millimetre emission lines obtained with EMIR at 
the IRAM 30-m telescope. Their dust and gas properties provide 
firm evidence that they are indeed gravitationally lensed galax¬ 
ies, as is further supported through interferometric observations 
of their dust and gas morphologies, already obtained for most 
sources. 

We used the WISE survey at 22 /am and the 1.4 GHz VLA 
FIRST survey to show that the far-infrared continuum of the 
Dusty GEMS is not dominated by the radiation of powerful 
AGN. In particular, we find that buried quasars cannot make a 
dominant contribution to their observed FIR SEDs. All SEDs 
are well fitted with a single modified blackbody distribution 
with temperatures T& = 33-50 K, covering the range of high- 
redshift starburst galaxies, as well as more gradually, but still 
intensely star-forming, high-z galaxies on the “main sequence.” 
They show a wide scatter about the local far-infrared radio cor¬ 
relation, with ^-parameters ranging from 2.0, as has previously 
been found for high-z galaxies by some authors, to above 2.7, 
which suggests a considerable excess of FIR relative to syn¬ 
chrotron emission. One plausible interpretation is that this is 
probably a consequence of their turbulent ISM, but this needs 
to be confirmed through more detailed studies comparing the re¬ 
solved radio emission with other source properties. All galaxies 
have gas-to-dust ratios of 40 - 140, consistent with a low CO- 
to-H 2 conversion factor, as expected for massive, dusty starburst 
galaxies with metallicities above solar. A full analysis of the spa¬ 
tially resolved properties of these galaxies, as well as detailed 
lens modeling, is on-going. 

Strongly lensed high-redshift submm galaxies represent an 
excellent opportunity to study gas heating and acceleration, and 


the mechanism driving star formation in the most vigorous star- 
bursts in the early Universe. Detailed observations of the dust, 
stellar populations, and multiple emission and absorption lines, 
in particular with submm and millimetre interferometry have al¬ 
ready been obtained and will be discussed in future papers. 
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